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ABSTRACT 

Formulas  are  deduced  for  vibrating  systems  of  one,  two, and 
three  dimensions.  Undamped  and  damped  free  vibrations  and 
harmonic  forced  vibrations  are  treated.  Methods  are  proposed 
for  calculating  the  damping  constants  from  test  observations 

I.  INTRODUCTION 

In  developing  formulas  for  vibration  and  possible  flutter  of  structures 
1/2 

such  as  rudders,  '  it  may  be  necessary  to  include  damping  forces.  Since 

these  forces  are  not  easy  to  calculate,  methods  of  determining  them  from  test 

3* 

observations  may  be  needed.  The  basic  theory  for  two-  and  three- 
dimensional  cases  will  be  considered  and  feasible  methods  of  observation 
will  be  sought.  First,  however,  formulas  for  the  one-dimensional  system 
will  be  written  to  assist  in  treating  the  main  problem.  For  convenience  of 
reference,  a  summary  of  the  results  is  given  in  Table  1;  see  pages  30  and  31. 
II.  ONE-DIMENSIONAL  VIBRATIONS 
Assume  as  the  equation  of  motion 

mx  +  cx  +  kx  =  P(t)  W 

in  which  m,  c,  and  k  are  positive  constants,  x  =  dx/dt,  and  P(t)  denotes  an 
applied  force  varying  with  the  time  t. 


1  References  are  listed  on  page  33. 

*  In  Reference  1  (see  pages  78  and  83),  certain  damping  terms  were  omitted 
from  the  flutter  equations  because  methods  for  determining  these  terms  from 
experiments, were  unknown  to  the  authors  at  that  time.  These  flutter  equa¬ 
tions  including  the  damping  terms  originally  omitted  are  of  the  same  form  as 
the  equations  given  here  for  the  three-dimensional  case. 
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1 .  DAMPED  FREE  VIBRATIONS 


If  P  =  0 ,  the  general  solution  of  Equation  [l]  can  be  written  (as  is 
easily  verified)  as  follows  in  terras  of  independently  arbitrary  amplitudes  a 


and  b: 


2  —  ■£ 

If  c  <4  mk  (less  than  critical  damping):  x  =  e  ^  (a  cos  cot  +  b  sin  cot) 


r  9  V  1  9  9 

where  and  to  =  -  -  -7  m  e 

n  2m  m  4 

If  =  4  mk  (critical  damping):  x  =  (a  +  bt)e  ^ ,  u =  — 

9  2m 

If  o'  >  4  mk  (greater  than  critical  damping):  x  =  ae  +  be 
where  jll^  and  ^  denote  the  following  two  values: 


^1,2  =2m  <C±^  C  ~4mk  > 


2  *  HARMONIC  FORCED  VIBRATIONS 

With  P  =  p  cos  <o  t  in  terms  of  arbitrary  constants  p  and  <o 

x  =  a  cos  to  t  +  b  sin  to  t 
o  o 


(k  -  mco  2)  +  c2co  21 

L  o  o  J 

(k  -  mw  2)2  +  c20.  2] 
L  o  o  J 


a  =  (k  -  mco  )  p 
o  '  K 


b  =  cco  p 
o 


2^  2  2  /  2  ,2*  2  a 

(k  -  mco^  '+ca)0  (a  +  b  )  ■  p  ,  ^ 


k  -  mco 


Thus  a  =  0  and  the  vibration  is  in  time  quadrature  relative  to  P 

when  co_=  k/m ,  which  is  the  value  of  co  for  undamped  free  vibration.  The 

J  2  2 

maximum  amplitude  or  maximum  of  \|  a  +  b  forgiven  p,  however,  occurs 


when 


(d/dco  )  (k  -  mco 


2  2  2  2] 

’o  >  +  °  “oJ  = 


or  when 


! 


2  k  c2 
w  =  — - r 

°  m  2m2 

2 

This  differs  from  k/rn  by  twice  as  much  as  does  to  in  a  damped  free 
oscillation. 

These  formulas  exhibit  several  features  for  which  analogs  mayreason- 
ably  be  expected  in  more  complicated  cases,  namely: 

(1)  Two  independent  modes  of  damped  free  vibration  occur.  Their 
amplitudes  can  be  chosen  to  make  x  and  x  agree  with  any  assumed  initial 
values. 

(2)  These  free  vibrations  are  oscillatory  provided  the  damping  constant 
c  is  not  too  large;  in  this  case,  c  produces  only  a  second-order  change  in  the 
oscillatory  frequency. 

(3)  In  a  harmonic  forced  vibration,  c  introduces  a  component  of  x  in 

time  quadrature  relative  to  the  applied  force  P  (proportional  to  sinco^t 

instead  of  to  cosw  t). 

o 

(4)  x  is  entirely  in  quadrature  relative^to  P  when  the  forcing  frequency 

factor  go  eauals  the  value  of  co  for  undamped  free  vibration, 
o 

(5)  The  maximum  amplitude  of  x  for  forcing  at  given  p,  when  damping 
is  present,  occurs  at  an  differing  from  the  undamped  free  <x>  by  more  than 
does  the  oscillatory  w  in  damped  free  vibration. 

3 .  EXPERIMENTAL  DETERMINATION  OF  c 

If^r^zO,  its  value  can  easily  be  determined  from  a  curve  showing 
either  x  or  x  as  a  function  of  t  during  damped  free  motion.  Then  c  =  2m yj. . 
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■JW  ... 


If  O)  is  also  determined  from  the  curve,  the  ratio  k/m  can  be  calculated  as 
2  2 

k/m  =  co  +A4  .  To  determine  k  and  m  separately,  one  of  them  must  be 
known  from  some  other  source. 

Or,  during  a  damped  forced  vibration  the  ratio  b/a  may  be  observed  as 
the  ratio  of  the  components  of  x  respectively  in  lagging  quadrature  to  P  and 

o 

in  phase  with  p,  or  the  equal  ratio  for  x.  (Note  that  here  x  =  x).  Then 

1  2*  b 

c  -  —  (k  -  mw  )  — 
co  o  a 

o 

In  this  case,  the  values  of  both  k  and  m  must  be  known. 

III.  TWO-DIMENSIONAL  VIBRATIONS 

Assume  that  the  kinetic  energy  T  and  potential  energy  V  of  a  two- 

4  S* 

dimensional  system  can  be  written  as 

T=iml*2+  2  m2/+  ml2^'  v=2  klx2+2k2y2  +  k12Xy 

in  which  x  and  y  are  generalized  coordinates  and  m.  ,  m  ,  m  are  inertial 

A  6  X  6 

and  k^,  k^,  k^  elastic  constants,  of  which  only  m^  and  k  £  may  be 
negative.  Substitution  of  first  q  =  x  and  then  q  =  y  in  Lagrange's  equation 
or 

d  aT  ,  av  -0 

dt  aq  aq 

gives  as  equations  of  motion 

ml^i  +  klX  +  ml2^  +  k12y  =  ^ 

ml2Vk12X+m2y  *  V  =^{t) 

.in  which  P(t)  and  Q(t)  represent  the  total  generalized  forces  acting  on  the 
*  Also  see  Appendix  A  of  Reference  1 .  * 
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system  (not  including  internal  elastic  forces).  Part  of  P  and  Q  may  be  due 
to  linear  damping  forces.  Expressing  the  latter  in  terms  of  damping 

constants  c^,  c^t  c^'  C2i'  * 'aQ  equations  of  motion  may  be  written: 

♦ 

+  klX  +  mi2V  +  k12y  +  V  +  °12y  =  P(t)  [2] 

mi2^  +  k12X  +  m2^  +  k2y  +  C21X  *  c2y  = 
in  which  P  and  Q  represent  possible  external  forces  acting  on  the  system 

(aside  from  damping  forces). 

Certain  restrictions  on  the  possible  values  of  the  constants  are  worth 
noting.  Let  x  and  y  be  so  chosen  that  T  and  V  are  never  negative.  Damping 
effects  can  never  increase  the  sum  T  +  V.  Multiply  the  first  of  Equations  [2] 
by  x  and  the  second  by  y  and  add  the  two  equations.  The  sum  of  the 
resulting  m  and  k  terms  is  easily  seen  to  equal  (d/dt)  (T  +  V);  hence,  if 


P  =  Q  =  0 


dT  (T  +  v)  =  -  cl*2  -  c2y2  -  (C12'+  c2l]  *y 


To  keep  (d/dt)  (T  +  V)  from  ever  being  positive,  it  is  necessary  that  >  0, 
c^  >  0,  since  either  x  or  y  may  vanish.-  Similarly,  to  keep  T  >  0  and  V  >0, 
it  is  necessary  that  m^,  m^,  k^,  and  k^  all  be  >  0. 

Further  restrictions  may  be  inferred  from  the  following  theorem.  LetoC, 

(3,y ,  e,  g  be  real  numbers.  Then 

2  2  2  2 

+  /3g  .+  /eg >0  or  de  +  J3g  >-ye g  fgi 


for  all  values  of  e  and  g  if  and  only  if 


<^>0,  /3>0,  y  <  8 


To  prove  this,  note  first  thatch  and  j8  cannot  be  negative  because  of 


cases  in  which  only  e  =  0  or  g  =  0.  Relation  [3]  then  clearly  holds  if  e  and 
g  are  such  that  /eg  >  0. 

Suppose,  however,  that  ye g  <  CL  Then  Equation  [3]  in  its  second 
form  is  equivalent  to  the  following: 

2 

(cCe2  +  fig2)  >  (ye g)  [3a] 

provided  that  positive  square  roots  are  taken  in  passing  back  from  Equation 
[3a]  to  Equation  [3]  .  But 

{c(e2+f3g2)  =  Ke2-/3g2)  +  4*/0(eg)2 

2  2 

Hence,  ifcOO  and/3>0  and  if  e  and  g  are  chosen  so  that  °^e  =  fig  » 

2  p  2  2 

then  -  4^/3  (eg)  .  Thus  Equation  [3a]  can  hold  generally  only  if 

either  or  yS  vanishes ,  Equation  [3]  requires  that  y-0 .  Conversely 

2  2  2  2  2. 2 
if  the  condition  that  4 r/.fi^ry  is  met  butene  =£  /3g  ,  then  (c4e  +  fig  ) 

>4cCfi(e g)2  >y2(eg)2  and  Equation  [3a]  holds,  also  Equation  [3]  * 

Substitute  here^rm^/2,  fi-m^/2,  y-m^,  e  =  x  and  g  =  y;  next, 

<=>£  =  k^/2 ,  fi-k^/2,  y- k12'  e  =  x  and  g  =  y;  and  finally  ^ ,  fi- c^, 
y  =  c  +  c  ,  e  =  x  and  g  =  y.  Compare  the  resulting  expressions  with 
expressions  previously  written  for  T,  V,  and  (d/dt)  (T  *  V),  It  will  then  be 
clear  that,  to  prevent  T  and  V  from  ever  becoming  negative  or  (d/dt)  (T  +  V) 
positive,  it  is  necessary  and  sufficient  that 

ml2  “  mim2'  k12  -  klk2'  ^C12  +  C21^  -  4ciC2  M 
These  restrictions  will  be  assumed  to  hold. 

It  follows  then  also  that 


2m12k12  ^  mlk2  +  m2kl*  C12°21  -  C1C2 


[5a.  b] 
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~«.A< 


For  (mjkg  ♦  mgkj)2  =  (mjkg  -  m^)2  +  ^m^k^  >  ^^k^ 

by  relations  [4]  .  (Note  that  a  square  cannot  be  negative.)  Similarly,  in 
any  case  4c12c21  <  4c12c21  +  (c^-c^)2;  (c^  +  c2i)2;  hence,  by 
Equation  [4]  ,  4c12c21<  40^2* 

Two  other  relations  that  can  be  inferred  in  a  similar  way  from  relations 


(C12+  C21)ml2<Cim2+  W  (C12  +  C21)k12<  Clk2  +  C2kl 
1 .  UNDAMPED  FREE  VIBRATIONS 

Undamped  free  oscillations  merits  consideration  as  background  for 

study  of  the  damped  case.  Let  c  =  c  =c  =  c  =0,P=Q  =  0.  Then 

I  z  1  z  z  1 

Equations  [2]  become 

n^x  +  k^x  +  m  y  +  k^y  =  0  m^x  +  k^x  +  n^y  +  k2y  =  0  [6] 

Two  special  cases  may  first  be  noted.  According  to  Equations  [6]  ,  x 

2  2 

can  vibrate  while  y  =  0  only  if  k^  -  m^co  and  k^2  -  m^co  are  k°*h  zero „ 
The  first  condition  fixes  co  at  ^jk^/m^;  the  second  requires  that  either 

ml2  =  ^12  =  ^  or  mi^i2  =  mi2^1‘  Similarly,  y  can  vibrate  with  x  =  0,  and 
co  =^[  ^2/m2  only  if  either  m12  =  k  2  =  0  or  m2k12  =  m12k2* 

If  x  and  y  vibrate  together  in  proportion  to  coscot,  the  following 


equations  must  be  satisfied: 

2  2 
(kx  -  m1co  )  x  +  (k12  -  m12co  )  y=  0 

2  2 

(k12  ~  mi2w  )  x  +  (k2  "  m2C0  )  7  =  0 

Elimination  of  x  and  y  gives  for  the  determination  of  co  the  following  equation 

2  2  2  2 
(k1  -  m1(j  )  (k2  -  m2u  )  -  (k12  -  m^co  )  =  0  [7a] 
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I  : 


.1  l 

)  \ 


if  l 

%  \ 


1  ♦ 

j '.}  ( 


-  _ 


or 


(mim2  ~  mi22^4  “  ^mlk2  +  m2kl  ~  2mi2ki2^2  +  ^klk2  ”  ^12  ^  =  0  t7b] 

2  2 
Ifk12  =  kIk2'  one  root  of  Equations  [7  b]  is:  co  =  0.  Alternatively, 

o 

if  m^m2  =  mi2  •  only  one  mode  of  vibration  is  possible. 

2  2  2 
Assume  now  that  <  k^k^  and  m^  <  m^m^.  To  ^oca^e  w  / 

consider  L,  the  left-hand  member  of  Equation  [7a]  or  Equation  [7b]  ,  as  a 

2  2  2  2 
function  of  co  .  At  to  =  0#  L  >  0;  but  when  co  has  increased  to  co 

min 

representing  the  lesser  of  the  two  values  k^/rn^  and  k^/m^,  then  it  is  clear 

2 

from  Equation  [7a]  that  L  <  0.  Hence  L=  0  at  some  positive  value  of  co 
2 

less  than  to  .  .  Also  at  the  greater  of  the  values  k  /m  and  k  /m  ,  L<0, 
mm  —  11 

2  ' 
but  as  to  co  it  is  clear  from  Equations  [7b]  that  L>  0.  Hence  a  second 

o 

root  of  Equation  [7a,  b]  occurs  at  a  value  of  to"  greater  than  both  k  /m  and 

VS- 

Thus  two  different  modes  of  vibration  of  the  system  are  possible  with 
both  x  and  y  vibrating .  •  In  each  mode 


Z  = 
x 


k12  -  ml2C°‘ 
,  2 
k2  ~  m2W 


k^  -  m^co 
k12  "  mi2W‘ 


2 .  DAMPED  FREE  VIBRATIONS 

Let  P  =  Q  =  0  so  that  Equations  [2]  read 

ml*+  klX+  mi2V  +  k12y  +  C1X  +  Cl2y  =  ° 

ml2X+k12X+  m2V  +  V+  C21X  +  C2y  =° 

In  special  cases  especially  if  m ^  =  k^  =  0  and  c^  =  0  so  that 
Equation  [8a]  reduces  to  Equation  [l]  with  P  =  0,  x  can  vary  while  y  =  0;  or, 


[8a] 

[8b] 
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similarly,  if  -  0,  y  alone  may  vary.  Such  cases  will  not  be  discussed 
further  here. 

For  the  general  case,  solutions  may  be  sought  in  which 

x  =  ae^,  y  =  be 

where  a,  b,  and  X  are  non-zero  constants,  real  or  complex.  Substituting  in 
Equations  [8a,b]  and  canceling  out  e 


(nij  X2  + 

kj  +  cAX  )  a  +(m12X  +k12  +  c12X)  b=  0 

[9a] 

(m  1 2  X  2+ 

k12  +  °21  X)  a  +  'm2  +  k2  +  °2  X )  b  =  0 

[9b] 

The  result  of  eliminating  a  and  b  from  these  equations  may  be  written: 

£4  X  4  +  e3  \3  +  e2  X2  +  e^X  +  eQ  -  0  [10] 

where 

eo  =  klk2  -  k122-  «1  =  Clk2  *  °2kl  -  (°12  *  °21)l!12 

*2  =Blk2*”2kl--Jml2k12*cl°2  -°12c21 

*3  =  °lm2  *  °2ml  '  (o12  *  c21)m12'  e4  =  ”lm2  -  "122 

The  coefficients  •  .  .  .  are  all  >  0,  according  to  Equations  [4] 

and  [5a,b,c,d]  .  Hence  no  root  X  of  Equation  [lO]  can  be  a  positive  real 

number.  Probably  if  the  damping  is  strong  enough,  negative  real  roots  may 

occur,  possibly  even  four  in  number,  but  this  difficult  question  is  of  little 

practical  interest  here. 

For  the  general  case,  write \=  -  fj.  +  iw  where  i  =  in  and  ^  and  co 
are  real  numbers.  The  following  two  equations  result  from  substituting  in 
Equation  [10].  then  equating  the  real  and  imaginary  parts  separately  to  zero; 
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A. 


and  dividing  the  imaginary  equation  by  iw  on  the  assumption  that  cj  z^z  0: 

4  2  2  2  3  4 

V*  “  ^2~3  e3^+  6  +  eo  "  ei/x+-2^  -<Vt  +  Vi  =  0  [lla] 

«j  -  Vo2  -  (2^  -  4  oj2 ) ^.  +  3e3/i2  -  4  e^/jL  =  0  [lib] 

2 

These  equations  determine^  and  co  .  The  conjugate  quantity  -  f_i  -ico  • 
is  then  also  a  root  of  Equation  [10].  Since  there  are  only  four  roots  in  all, 
there  can  be  only  two  pairs  of  values,^  and  to..  ,  and /.t  and  to  *  These 
pairs  define  two  modes  of  damped  oscillation.  Since  damping  cannot 
increase  the  total  eneigy,  it  must  turn  out  that  both  p  and  jj,^  are  positive. 

To  obtain  real  expressions  either  the  real  parts  of  all  quantities  (i.e., 
solutions)  may  be  chosen  or  the  imaginary  parts  divided  by  i;  the  two  pairs 
of  real  solutions  thus  obtained  are  in  relative  time  quadrature.  The  value  of 
the  ratio  b/a  for  each  mode  may  be  obtained  from  Equations  [9a, b]  .  Since 
usually  b/a  will  turn  out  complex,  there  will  generally  be  a  difference  of 
phase  between  x  and  y  as  functions  of  the  time. 

Thus  four  real  expressions  are  obtained  representing  four  independent 
damped  oscillations.  For  these  oscillations,  x  and  y  can  be  written  thus: 
x  =  e”^1*  (A1  cosc^t  +  A^  sino^t),  y  =  r^  e  ^^[a^  cos  (c^t  +  sin^t  +  4^)] 

or 

y  =  e  yU'2t  (A2  cosWjt  +  A^  sinc^t),  y  =  r2  e"^2*  [a2  cos  (o^t  +€  )+A‘  sin(«  t+e  )] 
Here  ,  A^ ,  A^  are  independent  arbitrary  constants  which  can  be 
adjusted  to  fit  any  assumed  initial  values  of  x,  x,  y,  y.  It  should  be  noted 
that 

10 


f 


i 


j 


•t 


/  2  2.  ,  ..  ,  2  2n 
X  =  “  (co1  X  and  y  =  -  fa  -/x^y 

in  any  one  mode  whereas  in  the  other 

"  /  2  2»  *»  <  2  2« 
x  =  -  (w2  "  X  dnd  7  r  "  ^2  2  7 

Only  small  damping  effects  appear  to  be  important  in  practice.  Hence 
no  general  discussion  of  Equations  [lla,b]  will  be  undertaken  here. 

If  the  c's  are  sufficiently  small,  ji  will  also  be  small,  and  the  co¬ 
efficients  and  are  likewise  small.  Consequently  all  terms  in 
Equation  [lla]  containing  are  small  at  least  to  the  second  order,  and  the  last 
three  terms  in  Equation  [lib]  are  small  to  the  third  order.  For  an  approximate 

solution,  these  terms  may  all  be  dropped.  Then  Equation  [lla]  becomes: 

4  2  r  n 

€4u  “  +  €o  =  ^*  This  agrees  with  Equations  [7a, bj  for  the  case  of 

no  damping  so  that  to  the  degree  of  approximation  under  discussion,  the 

oscillation  frequencies  are  the  same  as  if  there  were  no  damping.  From 

Equation  [lib]  the  approximate  value  of  /x  is 


ei‘V° 


V 2  V 


[12] 


More  accurate  solutions  can  be  obtained  from  Equations  [lla,b]  by  a 
process  of  successive  approximation. 

3 .  HARMONIC  FORCED  VIBRATIONS 

If  the  applied  forces  are  harmonic  functions  of  the  time  t,  they  cause 
harmonic  vibrations  of  x  and  y.  At  the  start  there  may  also  exist  superposed 
damped  free  oscillations  whose  amplitudes  can  be  adjusted  so  as  to  produce 


'Jf 


on  the  whole  any  initial  values  of  x,  y,  x,  y.  These  damped  free  oscillations 
will  be  assumed  to  have  died  out. 

Since  in  the  one-dimensional  case,  the  presence  of  damping  intro¬ 
duces  a  phase  difference,  assume: 


P  =  p  cos  co  t  +  p’  sin co  t, 
o  o 

■x  =  an  cos  to  t  +  a'  sinco  t, 
1  p  1  o 


O-  q  cos  co  t  +  q’  sinco  t 

°  [13] 

y  =  a„  cos  co  t  +  a'  sin  co  t 
1  2  o  2  o 


In  Equations  T?.!  the  cosco' t  and  sinco  t  terms  must  balance  separately. 

L*  J  '  o  o 

After  canceling  the  time  factors,  the  result  is  the  following  four  equations 
<kl  -  ml“o2)al  +  ClVi  +  {k12  -  ml2tJo2)a2  +-C12“oa2  =  P 


-°lVl  +  (kl  -  ml“c2)ai  "  C12woa2  +  (k12  ~  ml2Wo2)a2  =  P' 

(k12  "  ml2wo2)al  +  C2lVl  t  (k2  "  m2Wo')a2  +  c2Woa2  =  q 
‘C2lVl  +  (k12  "  mi2a,o2)ai  "  C2woa2  +  (k2  " 


[14] 


Here  p,  p',  q,  q',  a  ,  a^,  a  ,  are  eight  real  numbers.  In  general,  any 

four  of  them  can  be  assigned  arbitrarily;  the  equations  then  fix  the  values 

of  the  other  four.  Furthermore,  since  cosco  t  and  sinco  t  differ  only  in 

o  o 

phase,  the  zero  for  t  can  be  so  adjusted  that  any  chosen  one  of  the  eight 


quantities  a^  .  .  .  .  q'  vanishes,  without  altering  the  physical  form  of  the 
vibration.  Thus  all  cases  can  be  covered  while  keeping  one  coefficient 
zero. 

In  particular,  Equations  [l4]  may  be  solved  for  the  amplitudes  a^,  a^# 
a2'  a2  causec*  ^7  given  applied  forces  represented  by  p,  p’,  q,  q1.  The 
determinant  a  of  the  coefficients  of  a^,  a^#  a^,  a^  is  easily  found  to  have 
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the  value 


*{(kl  -  m1c°o2){k2  "  m2Wo2)  ~  (k12  "  ml2“o2)  ]  .  "  “o4(c1C2  “  C12°21)2 
If  there  is  no  damping,  comparison  with  Equation  [7a]  shows  that  a  =  0 


when  equals  the  value  of  co  for  either  of  the  frequencies  of  undamped  free 
vibration  of  the  system. 

If  ,  ci2 ,  c2i  are  merely  all  small,  a  will  vanish  at  two 

slightly  modified  frequencies  that  differ  also  slightly  from  the  frequencies 
of  damped  free  oscillation.  As  co^  approaches  either  of  these  frequencies  at 
which  a  —  0  while  p,  p' ,  q,  q*  remain  fixed,  the  amplitude  of  the  forced  vi¬ 
bration  becomes  large  (the  phenomenon  called  resonance). 

4.  EXPERIMENTAL  DETERMINATION  OF  c  ,  c  ,  c21 

One  method  is  to  make  "bumping"  observations  by  starting  a  motion 
and  recording  it  as  it  decays  .  By  proper  adjustment  of  the  initial  values  of 
x,  x,  y,  y,  the  system  can  be  made  to  vibrate  in  either  of  its  two  modes  of 

damped  free  vibration  with  the  other  mode  absent.  Observations  may  be 

2 

made  of  either  x  and  y  or  x  and  yas  functions  of  the  time  since  x=  -  (co^  - 

2  ..  2  2  ..  2  2  .. 

)  x  and  y  =  -  (u  -  /j,  )  y  in  one  mode  and  x  =  -  (a>2  -/^  x,  Y  ~ 

2  2 

-  (co2  -/x 2  )  y  in  the  other.  From  these  observations,  values  can  be  cal¬ 
culated  for  each  mode  of  the  frequency  u>,  the  damping  constant  fx  ,  and  the 
amplitude  ratio  r  and. phase  €  of  y  relative  to  x,  giving  the  eight  known 
quantities 


W1  “2  ^1^2  rl  r2  €1  €2 
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Insertion ~of  co  and/x,  and  then  of  to  and /x  for  to  and  /x  in  Equations 

XX  £*  t* 

[lla,b]  then  provides  four  equations  which  can  be  solved  numerically  for 

c  ,  c0 ,  c  ,  and  c0  provided  the  six  constants  m  ,  m  ,  m  ,  k  ,  k  ,  and 
1  Z  1Z  Z  x  1  Z  1Z  i  z 

k^^  are  known..  It  might  be  more  accurate,  however,  to  use  equations 
containing  the  constants  e.  and  e  which  differ  from  zero  only  because  of 

X  ^ 

damping.  If  bumping  observations  are  to  be  used,  further  study  of  the 
methods  of  calculation  should  be  made.  The  damping  may  be  weak  enough 
to  justify  the  use  of  simplifying  approximations. 

It  may  be  worth  noting  that  observation  of  all  eight  quantities  to 
should  make  possible  the  calculation  of  nine  of  the  ten  quantities  m^. 


l2'  mi2'  V  k2'  k12'  V  C2'  C1 2 ' 


and  c  •  F°r  a  restriction  exists  on 


the  possible  variation  of  these  quantities.  Let  Equations  [8a,  b]  be 
multiplied  by  an  arbitrary  constant  _s.  The  new  equations  may  then  be 
regarded  either  as  equations  in  a  different  form  for  the  original  system  or  as 
equations  for  a  different  system  having  constants  _s  times  as  great  but  the 
same  damping  modes  as  the  original  system.  In  order  to  know  which  system 

of  this  similitude  class  the  observed  constants  . refer  to, 

it  is  necessary  to  know  at  least  one  of  the  ten  quantities  . 

Then  the  remaining  nine  can  all  be  calculated  from  the  eight  observed 
constants  co^  .  .  .  .  .  .  .  . 

(If  Equations  [8a,  b]  are  multiplied  by  different  numbers,  they  are  still 
valid  for  the  original  system  but  cannot  be  regarded  as  equations  in  the  same 


form  as  Equations  [8a, b]  for  a  different  system  because  the  new  m12^  m21 
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and  k^^b  k21*^  Even  ^  the  initial  values  of  x>  x/  Y*  V  cannot  be  properly 

adjusted,  since  one  mode  will  usually  die  out  before  the  other,  both  sets  of 

values,  ^  and  to^  and  jx^  and  co^ ,  can  be  inferred  from  the  same  curve  of  x 

or  y  as  a  function  of  time.  If  both  modes  persist,  it  is  still  possible  to 

7-10 

observe  each  mode  in  turn  by  means  of  a  filter.  Or  a  vibrator  may  be 

used  and  adjusted  in  frequency  so  as  to  be  in  resonance  with  one  mode; 
then,  after  the  vibrator  is  removed,  a  damped  free  oscillation  will  occur  in 
this  mode  only. 

If  c^  -  °2\  ~  ci  anc*  C2  can  calculated  from  /x  ^  and fx 
Otherwise  the  observed  values  of  /a  ^  and  ^  furnish  only  two  relations 
among  the  four  quantities  c  ,  c^,  c  ,  c  . 

1  A  i  6  6  1 

An  alternative  method  is  to  study  forced  harmonic  vibrations  produced 
by  applied  forces  P  and  Q  whose  relative  amplitudes  and  phases  can  be 
controlled.  (Applied  forces  are  pure  P  when  they  do  no  work  during  varia¬ 
tion  of  y  alone,  or  pure  Q  when  no  work  is  done  during  variation  of  x  alone.) 
Two  alternative  procedures  will  be  described  which  require  no  measurements 
of  P  or  Q.  The  constants  ,  m2,  m^,  k^ ,  k2,  k  2#  however,  must  be 

known.  Either  x  and  y  or  x  and  y  may  be  observed  since  in  forced 

2  ..  2  2 

oscillations  x  =  -co  x,  y  =  -<o  y  and  co  will  be  seen  to  cancel  out  in  all 
o  o  o 

final  formulas. 

First  Procedure:  Isolation  of  c^,.  c2#  c  2,  c  in  turn.  Make 
observations  as  follows: 


(1)  Cause  x  to  vibrate  with  y  =  0.  Assume  pr  =  0,  so  that  a^  denotes 
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the  amplitude  of  the  component  of  x  that  is  in  phase  with  P  and  a[  the 

amplitude  of  the  quadrature  component  of  x.  To  do  this,  apply  P  =  o  cos  to  t 

o 

and  adjust  the  amplitude  and  phase  of  Q  so  that  a  =  a’  =  0.  Then  Equations 
[14]  reduce  to 


2 

(k_  -  m  )  a.  +  cw  a*  z  p 
1  1  o  1  1  o  1 


clVl  +  (kl  -  mlwo  )  ai  = 0 


(k12-mi2“o  )al  +c2lVl  =q 


The  magnitude  of  cjq  should  be  made  quite  different  from  ^  k^/m^ .  Only  the 
ratio  a^/a^  needs  to  be  observed. 

Probably  the  adjustment  of  Q  can  be  effected  most  conveniently  by 
varying  its  amplitude  (q)  +  (q‘)  until  a g  (or  the  component  of  y  in  phase 
with  P)  is  zero,  then  varying  the  phase  of  Q  (thus  varying  q‘)  until  the 
quadrature  amplitude  a^  of  y  equals  zero,  and  repeating  these  adjustments 
in  turn  until  both  a 2  and  a^  remain  negligibly  small. 

Then 


-  1  tv  2. 

C1  -  {kl  "  ml“o  } 

o  1 

(2)  Similarly,  to  keep  x  =  0,  apply  Q  =  qcos^t,  hence  q'  =  0,  and 
with  not  near  ,  adjust  p  and  p‘  so  that  a^  =  a^  =  0,  and  read 


a2^a2*  Then 


C2  =  Zir  (k2  -  m2Wo2) 

o  2 
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(3)  Cause  x  and  y  to  vibrate  in  phase  with  P;  that  is,  writing 

P  =  pcoscoot  with  p'  =  0,  adjust  q  and  q'  so  that  a'  =  a'  =  0.  Read  a  /a  . 

X  Cm  JL  Cm 

Then  from  the  second  one  of  Equations  [l4] 

al 

In  this  case  the  simplest  way  to  effect  the  required  adjustment  of  Q 
might  be  to  vary  its  amplitude  so  as  to  reduce  the  larger  of  a*  and  a*  until 

X  6 

a'l  =  d2#  t*ien  adiust  t^ie  Phase  of  Q  so  as  to  minimize  a*  and  repeat  these 

Cm 

adjustments  until  a^  and  a^  have  been  made  sufficiently  small. 

(4)  Cause  x  and  y  to  vibrate  in  phase  with  Q,  assuming  q1  =  0.  Adjust 
p  and  p'  so  that  aj  =  a^  =  0  nearly  enough.  Read  a^/a  .  Then  from  the 
fourth  of  Equations  [l4] 

c  -3C 
21  C2 

This  procedure  should  yield  the  most  accurate  values  of  the  four  c's, 

but  the  experimental  adjustments  required  may  be  considered  too  tedious. 

Second . Procedure :  Single- phase  forcing .  Apply  p  and  Q  in  any  known 

ratio  but  in  the  same  phase.  Write  P  =  p  cos  co  t,  Q  =  q  cos  to  t,  so  that 

o  o 

p*  =  q'  =  0.  Read  a^  a2  as  amplitudes  of  inphase  and  a'  ,  a*  as  amplitudes 

X  z 

of  quadrature  components  of  x  and  y.  Repeat  with  a  different  ratio  Q/P, 
distinguishing  the  amplitudes  thus  obtained  by  a  bar. 

Substitute  each  set  of  a's  in  turn  into  the  second  and  fourth  of 
Equations  [l4]  ,  in  which  p'  =  q'  =  0.  The  resulting  equations  can  be 
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•j 

t  i 


\ 

■  sy 


kA.... 


written: 

Wl  '  "oa2°12  =  (kl  '  ml“o2)  al  *  <k!2  -  m12"oZ)  a2 

VlC2  *  "oa2c!2  =  <kl  -  "W’ ai*  <k12  ‘ 

Vlc21*  “oa2c2  =  (k12  "m12“o2)  ai  *  <k2  -  m2“o2)  »2 
Vlc21  *  “0V2  =  (k12  -  m12"o2>  »1  *  <k2  -  m2”o2)  a2 

These  two  pairs  of  equations  are  easily  solved  for  c  ,  c  ,  and  c  ,  c 

1  12  2  21 

IV.  THREE-DIMENSIONAL  VIBRATIONS 
Let  x,  y,  z  denote  the  three  displacement  variables,  for  example  v;  y 
oC  motion  of  a  rudder  (see  Reference  1).  Then  linear  equations  of  motion 
can  be  written  as  follows: 

nyx  *  kjX  ♦  m12y  +  k^y  +  m^z  ♦  k^z  +  +  c^y  +  c^i  =  P(t)  [l5a] 

m12k*  +  k12X  +  m2*y  +  k2y  +  m23*Z*  +  k23Z  +  °21X  +  C2y  +  °23Z  =  QW  C15b] 

mi3k  +  k13X  +  m23y  +  k23y  +  m3Z  +  k3Z  +  °31X  +  °32y  +  C3Z  =  R(t)  C15cl 
Here  P,  Q,  and  R  are  generalized  external  forces  so  defined  that  the  rate  at 
which  they  do  work  on  the  system  is  always  Px  +  Qy  +  Rz.  The  m‘s  are  of 
the  nature  of  inertial  constants  and  the  k's  of  elastic  constants. 

Then  there  may  be,  as  in  Equations  [lSa,b,c]  ,  nine  linear  damping 


constants  V  V  Cg.  c^,  c^,  c  .  c  c  ,  c  The  six 


cross 


constants  c^,  etc.,  will  be  limited  in  relative  size,  as  in  the  two- 
dimensional  case,  since  the  damping  necessarily  tends  to  decrease  the  total 
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I 

...  .  -  ...  . .  ...  ....  .  ...  . 

energy  T  +  V;  they  are  likely  to  be  relatively  small  and  may  be  negligible, 
but  this  cannot  be  assumed  to  be  true  in  general  because  the  magnitudes  of 
all  nine  constants  will  vary  with  the  choice  of  the  variables  to  be  called  x, 
y,  z. 

The  situation  will  be  analogous  in  general  to  that  for  two  dimensions. 
If  P  -  Q-R-0  and  all  c's  are  zero,  there  will  be  solutions  of  Equations 
[l5a,b,c]  representing  three  modes  of  undamped  free  vibration.  If  any  c's 
do  not  vanish,  these  modes  become  three  modes  of  damped  free  oscillation; 
or,  if  the  c's  are  sufficiently  large,  one  or  more  modes  may  be  replaced  by 
two  modes  of  exponential  decrease  without  oscillation,  such  as  were 
represented  by  formulas  in  the  one-dimensional  case. 

In  the  oscillatory  case,  on  the  other  hand,  there  will  be  three  damping 
constants  /jl^,  In  any  one  mode  of  damped  oscillation,  the  three 

variables  x,  y,  and  z  may  be  assumed  to  be  proportional  to  e  ^1*  cos 
(co^t  +  <s),  in  another  mode  to  e  ^'2t  cos  (to^t  +6),  and  in  the  third  to  e  ^3* 
cos  (co^t  +  €),  the  phase  angle  €  being  different  in  general  for  x,  y  and  z 
and  different  in  the  three  modes. 

The  frequency  factors  will  not  be  quite  the  same  as  in  the 

undamped  vibrations,  but  the  difference  will  be  only  of  the  second  order  if 
the  damping  is  relatively  small. 

A  more  detailed  discussion  of  these  various  cases  follows; 
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1 .  UNDAMPED  FREE  VIBRATIONS 

If  P  =  Q  =  R  =  0  and  all  the  c*s  are  zero,  a  solution  of  Equations 
[l5a,b,c]  is 

x  =  a^  cos  tot,  y  =  a^  cos  cot,  and  z  =  a^  cos  cot,  a^,  a^/  and  a^  being  real 
numbers;  from  Equations  [l5a,b,c]  : 


(ki  "  miw2)  al  +  (k12  "  ml2w2)  a2  +  (k13  '  mi3C°2)  a3  =  ° 
(k12  “  m!2w2)  al  +  (k2  '  m2C°2)  V  (k23  “  m2 3 a3  =  ° 
(k13  ‘  ml3c°2)  al  +  (k23  “  m23°j2)  3 2  +  (k3  “  m3“?1  a3  "  ° 


Equating  the  determinant  of  a^,  a^,  a^  in  these  equations  to  zero  gives  the 
equation: 


-  mico2)(k2  -  m2co2)(k3  -  m^co2) 

+  2(k12  ‘  mi2tj2)(k13  "  ml3C°Z)(k23  ‘  m23‘°2) 

-  (kj  -  m1cZ)(k23  -  -  (k2  -  m2co2)(k13  -  m^J) 

2 

'  “  (k3  "  m3w2)(k12  "  mi2w2)  =°  C16] 


This  is  a  cubic  equation  in  to  whose  three  roots  furnish  the  frequencies  for 
three  modes  of  undamped  free  vibration.  Any  two  of  the  original  equations 
can  be  solved  for  the  ratios  of  a^ ,  and  a^  to  each  other  in  any  one  of 
the  three  modes  (see,  for  example.  Appendix  C  of  Reference  1). 


2 .  DAMPED  FREE  VIBRATIONS 

Assume  P  =  Q  =  R  =  0  and  write 


x 


Xt  Xt 

°ie  '  y  "  a2e 


Xt 

z  =a3e 
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where  a^,  and  X  may  all  be  complex  numbers.  Substitution  in 

Equations  [l5a,b,c]  then  gives: 

2  2 
X  +  c^X  )  +  (k^2  +  mi2  ^  +  ^12^^  ^2 


+  (k13+mi3X  +C13X,a3=° 


2  2 
^12  +  mi2  X  *  C21  X  )  ai  +  (k2  +  m2  X  +  C2X  ^  a' 


+  {k23+m23X  +C23X)a3  =° 


(k13  +  mi3  X  2  +  C31X  }  ai  +  (k23  +  m23  X  2  +  °32X  >  92 


(k3  +  m3  X  +  C3X)  a3  =  ° 

The  determinant  of  a^,  a2,  a3  in  these  three  equations  set  equal  to  zero 
gives  : 


2  2  2 
(k1  +  nij  X  +  cx  X  ) (k2  +  m2  X  +  c2X  ) (k3  +  m3  X  +  c3  X  ) 


+  (k12+mi2X  +c12X)(k23+  m53  ^  +cooX)(k1o+  m,oX2+C01X) 


12  23  23 


23  13  “‘13  31' 


+  (k12+m!2X  +C21X,(k23+  m23X  +  C32  X)  (k13  +  mi3  X  2  +  Ci3X) 

■(k1*n,lx2  *°lX)(k23*m23X2  *c23X)(t23‘  m23X  2  'C32X> 

(k2  ,m2X2  .c 2X)(k13.  n,13X2  *<=,3X)fc13*  ">i3*  2  <■  Cj.X) 


2  o 

~  (k3  +  m3X  +C3X  ^k12  +  mi2X  +cl2X^k12  +  ml2X  2  +  X  ^  =  0 


21 


[17] 


21 


) 


This  is  an  equation  of  the  sixth  degree  in  X  .  It  may  have  real  roots 
if  the  c's  are  large  enough,  perhaps  as  many  as  six  real  roots  .  On  the 
other  hand,  analogy  with  the  two-dimensional  case  suggests  that  if  the  c's 
are  not  too  large,  there  will  be  six  complex  roots  in  three  pairs:  --u^  +  ico^, 

"^2  iicV  3  iicV 

Two  equations  for  the  determination  of  and /x^,/x  /x 

analogous  to  Equations  [lla,b]  in  two  dimensions,  can  be  obtained  by  sub¬ 
stituting  X  =  ~}x  +  ho  and  separating  real  and  imaginary  parts.  In  the  three- 
dimensional  case,  however,  these  equations  are  voluminous  and  the  chance 
of  their  ever  being  put  to  practical  use  scorns  to  be  very  small,  hence  they 
will  not  be  written  out  here  in  full. 

For  practical  use  when  the  c's  and  hence  also  the  jx  *s  are  small, 
abbreviated  approximate  equations  can  be  obtained  by  omitting  all  terms  of 
second  or  higher  order,  that  is,  all  terms  containing  a  pov/er  of  jx  higher 

than  the  first  or  both  /x  and  one  of  the  c's  or  the  product  of  two  c's.  This 

2  2 

rule  of  approximation  justifies  replacing  X  in  Equation  [l 7\  by  -<*>  -2ioo/i. 

and  also  X  by  i<o.  Furthermore,  all  products  of  c  terms. may  be  omitted.  The 

first  of  the  six  products  in  Equation  [l7]  ,  for  example,  is  to  be  replaced  by 

2  2 
(k.^  -  nijW  -  2icom1  /x  +  icoCjH^  -  m^w  -  2iwm2  /x  +  iwcp 

2 

(IC3  -  m^co  -  2itom3  }x  +  i^c^) 

and  then  expanded  as 
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Z  2  2 

(j<^  ~  rn^co  )  0^2  *”  ^2^  )0^3  ~  m^co  ) 

+  iw  (-2m1  ft  +ci)(k2  -  m2cj2)(k3  -  n^oj2) 

2  2 

+  iw  (“2m2  /.l  +c2)(kf  -  m^co'Hkg  -  m3o  ) 

2  2 

+  ioj  (-2m3  /j.  -f-c3)(i<1  -  m^w  ){k2  -  m2<u  ) 

It  is  easily  seen  that  the  real  pari:  of  Equation  [17]  as  thus  reduced  is 
the  same  as  Equation  [16]  for  undamped  vibration.  Hence  the  frequencies 
of  oscillation  in  the  three  damped  modes  are  approximated  here  by  the 
frequencies  of  undamped  vibration  and  may  be  calculated  from  Equation  [16] 
To  shorten  the  notation,  write  now 

2  2  2 
Gj  -  k^  -  ,  G2  =  k2  -  m?w  ,  G3  =  k3  -  m3o> 

G12=  k12  _  mi2,c°  '  G13=  k13  "  m!3tj2'  G23=  k23  "  m23w2 


Then  it  will  be  found  that  the  imaqi: 


can  be  written  in  its  approximated  form  thus. 


Equation  [l  7]  divided  by  iw 


-2  ft  [m^GgGg  -  G232)  *  m2(G1G3  -  G^2)*  -  G^2) 


+  2ml2(G13G23  ~  G3G12)  +  2ml3  (G12G13  "  G2G13) 

+  2m23(G12G13  "  G1G230  +C1  (G2G3  "  G23  )  +  C2(G1G3  "  G13^ 
+  C3^G1G2  ”  G12  ^  +  ^C12  +  C21^G13G23  "  G3G12*  +  ^C13  +  C31* 


{G12G23  "  G2G13}  +  (C23  +  C32)(G12G13  “  G1G23}  =  0 
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*v 


2 

After  inserting  in  the  G’s  the  proper  value  of  to  for  any  one  of  the 
damped  modes,  this  equation  .is  easily  solved  for  an  approximate  value  of 
the  damping  constant  /.>,  for  that  mode . 

3  .  HARMONIC  FORCED  VIBRATIONS 


Assume 

P  =  p  cos  to  t  +  p*  sin  co  t,  Q  =  q  cos  co  t  +  q'  sin  to  t 

o  o  o  o 

R  =  r  cos  oj  t  +  r*  sin  co  t 
o  o 

where  p,  p1,  q ,  q’,  c,  r'  are  any  six  real  ainplitudes  and  co  is  any  positive 

real  number.  For  the  resulting  steady  vibration  write 

x  =  a  cos  w  t  +  a*  sinco  t,  y  =  a0  cosoi  t  •}  sin  co  t 

1010  2  o  2  o 

z  -  a0  cos  co  t  +  a*  sin  w  t 
J  o  3  o 

ai  *  *  *  •  a3  six  real  numbers. 


In  any  particular  motion,  by  a  proper  choice  of  the  origin  for  t,  any 


chosen  one  of  the  six  variables  P,  Q,  R,  x,  y,  z  can  be  made  to  vibrate  in 
proportion  to  cosco^t,  or  to  sinco^t.  Thus  any  one  of  the  twelve  amplitudes 
P/  P  ....  can  be  assumed  to  be  zero  without  altering  the  motion 

that  is  represented. 

Substitution  in  Equations  [l5a,b,cj  and  separation  of  sine  and  cosine 
terms  gives  six  equations.  To  shorten  the  notation,  write: 


F  =  k_  -  m_ cj 
1  1  1  o 


F2  ”  k2  m2a>o 


V  k3  "  'Vo' 


F12=  k12  ml2C°o 


F13_  k13  ml3wo  F23=  k23  "  m23Wo2 
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Then  the  six  equations  read: 

Vl  *  Wl  *  F12a2  *  °12"oa2  *  F13a3  *  °13“oa3  =  P 


[19a] 


■SVl  *  Flai  -  C12"oa2  "  ri2°2  -  °13“oa3  *  F13a3  =  P'  [l9b] 


F12al  *  °2l“oal  *F2a2  *  °2"oa2  ’  F23a3  *  C23“oa3  °  <* 

-21  Vl  *  F 1 2a  1  -  VoV  Vi  -  °23woa3  *  Va3  =  q'  ^ 

F13al  *  <=3lVi  ’  VV  C32«oa3 *  V»  *  Voa3  =  '  M 

-c3l"o“l  *  F13a'l  -  °32“oa2  *  F23a2  '  °3“oa3  *  Vi  =  r'  [l”] 
In  general  any  six  of  the  twelve  amplitudes  .....  r  can  be 
assigned  arbitrarily  and  the  equations  then  fix  the  values  of  the  other  six. 


4.  EXPESIMlMTAL  DETS'MIN  ATI  ON  OF 

C1.'>  c2}  c3>  c12’  C135  C2l’  C23’  °3l’  °32 

The  methods  described  for  a  two-dimensional  system  can  be  extended 
to  three  dimensions .  Determination  of  the  nine  damping  constants  from 
general  bumping  observations,  however,  will  not  be  discussed  here  because 
it  appears  to  involve  very  complicated  observations  and  calculations. 

A  feasible  alternative  might  be  to  lock  each  of  the  three  coordinates  in 
turn  so  as  to  hold  it  at  zero.  The  given  three-dimensional  system  could  thus 
be  studied  as  a  combination  of  three  two-dimensional  systems  and  the 
methods  already  described  for  such  systems  would  be  available. 

Of  three-dimensional  motions,  only  forced  harmonic  motions  will  be 


•ed  here  and  only  the  simplest  use  of  these,  J.n  such  motions,  x,  y, 

o  2  2 

o  ornin  1  s n g o'  i vg  1  v  to  ”0)  X/  —to  y*  —to  z  so  t licit  6ith.sr  x ,  Yt 

'  . '  '  4  o  o  o 


or  z  ma y  bo  rncasurod 


Apply  P ,  Q ,  11  *x>  any  ocnvon Ion'!;  ratio  but  all  in  tho  some  phase 
Assume  p’  =  q‘=  r!  =  0...  load  Lho  ro suiting  tnroe  inpbaco  arnpUtudoo  a.. 


a  and  tho  drone  cam 

o 


.mpHts.des  a.p  ,  the  latter  hdng 


relatively  enr.ll.  Repeat  with  different  ratios  of  P,  O,  R.  distinguishing 
the  a’s  thus  obtained  by  a  bar,  end  then  with  a  third  set  of  ratios,  mark! 
rl;e  a‘s  vr.-o  cl- >t >• ,  [\  possible  choice  would  be  to  use  cn.lv  P  the 
fhet  time,  only  O  o  s-.-co nd  K;rc-,  C',f  en.h,-  p  tine. 

Substitution  of  the  first  set  of  o berried  fa  ckuvblone  f 'i9b,d.?"l 
then  the  second  set  of  a's,  and  finally  the  thhj  sot  give?  tLrcn  groups  o 
three  equnuers  each  for  the  dote)  nd net  ion  of  '.ho  nine  c;s.  Stucco  in  all 


cases  p’ -  q‘ ~  r'  0,  the  equations  may  be  slightly  rearranged 
follows; 

W  a2°12  *  V'l3  =  «"  (Fl3i  +  FI2a2  +  FJ353) 

o 

a.c.  •>  a„c, +  a„c.  .  =  -----  (p  a'  +  F  s'  •;  F  a’ } 

1  i  2  H  3  13  o  v'  11.  12  2  13  3 

o 


T.  d  a 


c,  =  — -  fi-h a’  *  v  » 


6.  O  +  cl  C  <•  a  C  n  - --  fid  ^  -,  ■■  -,  t 

11  2  12  ’3  13  w  U  11  "  12-  1  ’  ' 

O 


t  Jib  "'3' 


a.(coi  +  a0c.  *  a  ,,c 


1^21  2^2  '  ~3~23  o)  v  12° I  '  ~  2~2  ‘  *  23’  o 

o 


Vh  "  V?:  V*s  '?  lFna‘l  -  F2a2  *r23ay 

o 

1 
} 

O 
1 


aiC2.1  *  a2C2  *  a3°23  r  coq  (F12ai  *  F232  +  F23a3* 

al°31  +  a2C32  +  a3C3  =  7T  <F.^  *  F-<  ai.) 
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CO,  13  1  "  23  2  3  3 


— i  A.- - 


aiC31  +  a2C32  +  a3°3  =  oj 


^F13al  +  F23a2  +  F3a3^ 


aiC31 *  a2C32  +  a3C3  =  ~  (F13aI  +  F23a2  +  W 

o 


Assuming  that  the  six  constants  Fj,  F^  F^,  F^,  F^,  F23  have  been 

calculated  from  the  constants  of  the  r.ys'em  and  the  chosen  value  of  co  ,  the 

o 


first  three  of  Equations  [20]  can  be  solved  for  c  ,  c  ,  c  ,  the  middle 

JL  X  &  lO 

three  for  c^,  c^,  c^^,  and  the  last  three  for  c0,  c^,  c^* 

The  computation  can  be  shortened  by  observing  differently.  Using 

chosen  p  and  q,  adjust  r  (p‘,  q',  r'  being  all  zero)  so  that  a,  =  0.  Read  a  , 

o  X 

a2'  ai/  a2‘  ar  Repeat  V7jth  a  different  pair  of  values  of  p  and  q,  making 

a3  =  0.  Read  a^,  a^,  a^,  a^ ,  a^.  Then  the  first  two  of  Equations  [20]  are 

easily  solved  for  c  and  c  9,  the  fourth  and  fifth  for  c  and  c  ,  and  the 
i  16  221 

seventh  and  eighth  for  c0.  and  . 

•  ox  32 

Repeat  using  two  pairs  of  values  of  p  and  r  and  adjusting  q  each  time 
so  that  a^  =  0.  Reada^,  a^,  a^,  a^,  a^  and  a^,  a^,  a^,  a a^.  Then 
the  same  three  pairs  out  of  Equations  [20j  yield  c  and  c  ,  c  and  c  , 

J.  lo  Z 1  2o 

and  and  c^ . 

All  of  the  c's  have  thus  been  obtained,  with  duplicate  values  of  c^ , 
c2i  and  c  .  Other  combinations  of  p,  q ,  r  may  be  used  in  a  similar  way. 

It  will  be  noted  that  neither  procedure  requires  actual  measurements 
of  P,  Q,  or  R. 


tne  same 
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I -  •  2 

to cc  l /  X  ;  /^cc  lA  ;  »#*>uoc  to 

Or,  if  c  does  n  o  f  change  when  X  =9~  1,  then  approximately  when  fx  is  small 

i  /  ^|T~  t  /  v3  .  6 

cj  cc  1/  X  ;  jj,  cc  1/  a  ;  .  ./xcc  co 

In  both  cases /.l  and  co  both  increase  if  X  <  1  and  decrease  if  X  >1,  thus 
varying  "in  the  same  direction." 

2  .  CONTRARY  MODES  FOR  A  GIVEN  SYSTEM 

Since  higher  frequency  tends  to  mean  higher  velocities  at  a  given 
amplitude,  it  might  reasonably  be  guessed  that  the  damping  will  be  greater 
in  modes  of  higher  frequency.  This  is  not  necessarily  the  case,  however, 
because  the  components  of  displacement  are  in  different  ratios  in  different 
modes  and  some  components  may  be  damped  more  heavily  than  others. 


As  a  simple  example,  suppose 

nrx  ,  kxx  +  k12y  ♦  Cji  =  0  m^y  +  k2y  -  k^x  ♦  =  0 

where  but  ^/m^. 

If  k,  =0  and  c.  =  c  =  0,  then  in  one  mode,  x  vibrates  with  y  =  0;  in  the 
other,  y  vibrates  at  much  higher  frequency  with  x  =  0.  If  kl2  =  0  but  c^ 
and  c2  are  merely  small,  then  the  two  frequencies  are  little  altered  by 
and  o  ,  and  the  damping  will  be  much  less  for  the  second  or  y  vibration 
than  for  the  first  or  x  vibration. 

Thus  higher  frequency  is  accompanied  here  by  lower  damping.  This 
conclusion  will  nc*  be  altered  if  k^2  is  merely  kept  small  but  not  zero,  so 
that  y  vibrates  a  little  in  the  first  mode  and  x  vibrates  a  little  in  the  second 
mode. 
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b  ration  occur  with  y  tnd  x  S 


jj<i>  HARMONIC  FORCED 


or  complex  (Xn«catlv*  ®r  complex) 


Table  1  -  Summary  of  Results  (Contin 
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equ&tton  In  X.  j  bratlo.-,*  cbs*in«d  from  t'-#  cab'.c 


For  each  mod#  or  value  of  w 
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